arXiv: 1508.03221 vl [cond-mat.str-el] 13 Aug 2015 


Permutation blocking path integral Monte Carlo approach to the uniform electron gas 

at finite temperature 

Tobias Dornheim, 1 ’* Tim Schoof, 1 Simon Groth, 1 Alexey Filinov, 1,2 and Michael Bonitz 1 

1 Institut fur Theoretische Physik und Astrophysik, 

Christian-Albrechts-Universitat, Leibnizstrasse 15, Kiel D-24098 , Germany 
2 Joint Institute for High Temperatures RAS, Izhorskaya Str. 13, 125412 Moscow, Russia 

(Dated: August 14, 2015) 

The uniform electron gas (UEG) at finite temperature is of high current interest due to its key 
relevance for many applications including dense plasmas and laser excited solids. In particular, 
density functional theory heavily relies on accurate thermodynamic data for the UEG. Until recently, 
the only existing first-principle results had been obtained for N = 33 electrons with restricted path 
integral Monte Carlo (RPIMC), for low to moderate density, r s = r/a b > 1. This data has been 
complemented by Configuration path integral Monte Carlo (CPIMC) simulations for r s < 1 that 
substantially deviate from RPIMC towards smaller r s and low temperature. In this work, we 
present results from an independent third method -the recently developed permutation blocking 
path integral Monte Carlo (PB-PIMC) approach [T. Dornheim et al, NJP 17, 073017 (2015)] which 
we extend to the UEG. Interestingly, PB-PIMC allows us to perform simulations over the entire 
density range down to half the Fermi temperature (9 = UbT/Ef = 0.5) and, therefore, to compare 
our results to both aforementioned methods. While we End excellent agreement with CPIMC, where 
results are available, we observe deviations from RPIMC that are beyond the statistical errors and 
increase with density. 

PACS numbers: 05.30.Fk, 71.10.Ca 


I. INTRODUCTION 

Over the last years, there has been an increasing inter¬ 
est in the thermodynamic properties of degenerate elec¬ 
trons in the quantum mechanical regime. Such infor¬ 
mation is vital for the description of highly compressed 
matter [1-3], including plasmas in laser fusion experi¬ 
ments [4-9] and in compact stars and planet cores [10- 
12]. In addition, the widespread density functional the¬ 
ory (DFT) approach crucially depends on the availabil¬ 
ity of accurate quantum Monte Carlo (QMC) data for 
the exchange correlation energy of the UEG, hitherto at 
zero temperature [13, 14]. However, in recent years more 
and more applications with highly excited electrons have 
emerged, which require to go beyond ground state DFT. 
Hence, there exists a high current need for an ab-initio 
thermodynamic description of the UEG at finite T. 

The widely used path integral Monte Carlo (PIMC) 
method, e.g. [18], is a powerful tool for the ab-initio sim¬ 
ulation of both distinguishable particles (often referred to 
as “boltzmannons”, e.g. [19, 20]) and bosons and allows 
for quasi exact results for up to N ~ 10 3 particles at finite 
temperature [21, 22]. However, the application of PIMC 
to fermions is hampered by the notorious fermion sign 
problem (FSP), e.g. [23], which might render even small 
systems unfeasible for state of the art QMC methods 
and is known to be NP-hard for a given representation 
[24] . With increasing degeneracy effects, permutation cy¬ 
cles with opposite signs nearly cancel each other and the 
statistical uncertainty grows exponentially. Hence, stan¬ 
dard PIMC cannot provide the desired results without 
further improvement. Brown et al. [25] have presented 
the first finite temperature results for the UEG down to 


r s = 1 using restricted PIMC (RPIMC) [26], a popular 
approach to extend PIMC to higher degeneracy, that is, 
lower temperature and higher density. To avoid the FSP, 
this method requires explicit knowledge of the nodal sur¬ 
face of the density matrix, which is, in general, unknown 
and one has to rely on approximations. The use of the 
ideal nodes for a nonideal system appears to be problem¬ 
atic, as has been shown for the case of hydrogen [27, 28]. 
In addition, it has been shown analytically that RPIMC 
does not reproduce the exact limit of the ideal Fermi gas 
(r s -A 0) [29, 30]. Therefore, the quality of the RPIMC 
data remains unclear. Indeed, recent configuration PIMC 
(CPIMC) [31, 32] results for the highly degenerate UEG 
by Schoof et al. [33] have revealed a significant disagree¬ 
ment between the two methods at small r s and low tem¬ 
perature. While the first application of a novel density 
matrix QMC (DMQMC) approach [34] to the UEG for 
four particles reports excellent agreement with CPIMC 

[35] , additional simulations of larger systems are needed 
to resolve the discrepancy towards RPIMC. For com¬ 
pleteness, we mention that QMC results by Filinov et al. 

[36] cannot be used as a benchmark due to the different 
treatment of the homogeneous positive background and a 
different account of the long-range Coulomb interaction 
[37, 38] than the usual Ewald summation. In this situa¬ 
tion an independent third first-principle method, capable 
to treat WDM parameters, would be highly desirable. 

In this work we, therefore, investigate the applicability 
of the recently developed permutation blocking PIMC 
(PB-PIMC) approach [39] to the uniform electron gas. 
The basic idea behind PB-PIMC is to combine antisym¬ 
metric imaginary time propagators [40-42], i.e., deter¬ 
minants, between all “time slices” with a higher order 
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II. THEORY 
A. Model Hamiltonian 

The uniform electron gas, often referred to as “Jel- 
lium”, is a model description of Coulomb interacting elec¬ 
trons with a neutralizing background of positive charges 
which are uncorrelated and homogeneously distributed. 
To describe an infinite system based on a finite number 
of particles, one implements periodic boundary condi¬ 
tions and includes the interaction of the N electrons in 
the main cell with all their images via Ewald summa¬ 
tion. Following the notation from [47], we express the 
Hamiltonian of the N electron UEG (in atomic units) as 


FIG. 1. Density-temperature plain around the warm dense 
matter (WDM) regime. PB-PIMC significantly extends the 
range of applicability of standard PIMC (qualitatively shown 
by the red dashe line) towards lower temperature and higher 
density while CPIMC is applicable to the highly degenerate 
and weakly nonideal UEG [33]. RPIMC data [25] are available 
for r s > 1. The orange area marks the conditions of WDM 
and inertial confinement fusion (IGF) [5]. 


factorization of the density matrix [43-46]. This means 
that each particle is represented by a “path” consisting 
of 3 x P coordinates (“beads”), where P is the number 
of high-temperature factors (or propagators). The appli¬ 
cation of determinants leads to a relieve of the FSP by 
an effective cancellation of positive and negative terms in 
the partition function, which belong to permutation cy¬ 
cles of different parity in standard PIMC. However, since 
the blocking is most effective if the thermal wavelength of 
a single propagator is of the same order as the mean in¬ 
terparticle distance, it is crucial to employ a higher order 
factorization scheme which allows for sufficient accuracy 
with only a few time slices. 
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Here R = niL and G = 112 /L denote the real and re¬ 
ciprocal space lattice vectors, respectively, with the box 
length L and volume V = L 3 . The specific choice of the 
Ewald parameter n does not influence the outcome of Eq. 
(1) and, therefore, can be used to optimize the conver¬ 
gence. PB-PIMC requires explicit knowledge of all forces 
in the system, and the force between the electrons i and 
j can be obtained from 

F ij = -Vi^( r i, r j) . (2) 

The evaluation of Eq. (2) is relatively straightforward 
and we find 


The details of our PB-PIMC scheme are described in 
section IIB, after a brief introduction of the employed 
model II A. In section III A, we present our simulation 
results starting with a detailed investigation of the con¬ 
vergence behavior with respect to the factorization of 
the density matrix. We proceed by simulating N = 33 
spin-polarized electrons, which is a commonly used model 
system of the UEG, see section IIIB. Interestingly, our 
PB-PIMC approach allows us to obtain accurate results 
over the entire density range and, therefore, to make a 
comparison with the pre-existing RPIMC and CPIMC 
results for the UEG. Finally, in section IIIC we investi¬ 
gate the applicability of our method with respect to the 
temperature. We find that PB-PIMC, in combination 
with CPIMC, allows for the simulation of the UEG over 
a broad parameter range, which includes the physically 
most interesting regime of warm dense matter, cf. Fig. 1. 
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B. Simulation method 

To calculate canonical expectation values with the PB- 
PIMC approach [39], we write the partition function in 
coordinate representation as 

z = E s S n ( CT ) / dR ( R l e ~ PH l^ R ) > ( 3 ) 

cr(=<Siv 
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with R = ri,...,rjv containing the coordinates of all 
electrons, Tt a denoting the exchange operator which cor¬ 
responds to a specific element a from the permutation 
group Sn and /3 = 1/k^T. For the next step, we make 
use of the usual group property of the density matrix 
in Eq. (3) and arrive at an expression for Z which re¬ 
quires the evaluation of P density matrices at P times 
higher temperature. However, instead of the primitive 
approximation e~ eH « e~ eK e ~ eV , with e = /3/P being 
the imaginary time step of a single propagator and the 
kinetic and potential contributions to the Hamiltonian I\ 
and V. respectively, we use the fourth order factorization 
[44, 45] 

e -tH r^, e ~v 1 eWa 1 £ -tieK e ~v 2 eW 1 -2a 1 

xe -t^K e ~vieW aie -2toek . ( 4 ) 

The W operators in Eq. (4) denote a modified potential, 
which combines V with double commutator terms of the 
form 

ilV : k},V} = _Y^\ Fl f , (5) 

TO z ' 


and, therefore, requires the evaluation of all forces on 
each particle, F. ( = — V*V r (R). Our final result for the 
partition function is given by 
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with the definition of the potential and force terms 
Vo = «iF(R a ) + u 2 y(R aA ) + v!V(R aB ) , (7) 
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with D being the dimensionality, see e.g. [40]. Eq. 
(6) contains two free coefficients, to and ai, which can 
be used for optimization, cf. Fig. 2, and the integra¬ 
tion is carried out over 3 P sets of coordinates, dX = 
dRo ... dRp_!dR 0 A ■ ■ ■ clRp—i^dRop . - • dRp_ip. In¬ 
stead of explicitly sampling each permutation individ¬ 
ually, as in standard PIMC, we combine configuration 
weights of both positive and negative sign in the de¬ 
terminants, which leads to a cancellation of terms and, 
therefore, an effective blocking of permutations. When 
the thermal wavelength of a single time slice, Xt ie = 
^2netih 2 /m, is comparable to the mean interparticle 
distance, the effect of the blocking is most pronounced 
and the average sign in our simulations is significantly 
increased. However, with an increasing number of prop¬ 
agators P, At 1£ decreases and, eventually, the blocking 


will have no effect and the sign converges towards the 
sign from standard PIMC. Hence, it is crucial to employ 
the high order factorization from Eq. (4), which allows for 
reasonable accuracy even for only two or three propaga¬ 
tors. We simulate the canonical probability distribution 
defined by Eq. (6) using the Metropolis algorithm [48] 
and refer to [39] for a more detailed description of the 
PB-PIMC method. 


C. Energy estimator 


The consideration of periodicity in the diffusion matri¬ 
ces requires minor modifications in the energy estimator 
presented in [39] , which can be derived from the partition 
function via the familiar relation 


1 dZ 
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Inserting the expression from Eq. (6) into (8) and per¬ 
forming a lengthy but straightforward calculation leads 
to 
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For completeness, we note that the total energy E splits 
into the kinetic and potential contribution, K and V, in 
precisely the same way as before [39]. 


III. RESULTS 
A. Convergence 

We begin the discussion of our simulation results by in¬ 
vestigating the convergence of the energy with the num¬ 
ber of imaginary time propagators P. To enhance the 
performance, the free parameters from the propagator, ai 
and to, can be optimized. In Fig. 2, we choose a\ = 0.33, 
which corresponds to equally weighted forces on all time 
slices, and plot the potential energy V, calculated with 
P = 2, P = 3 and P = 4, versus to over the entire 
possible range for a benchmark system of N = 4 spin- 
polarized electrons with 9 = 0.5 and r s = 4. To asses the 
accuracy, we compare these results with the exact energy 
known from CPIMC (green line). Evidently, the optimal 
choice for this free parameter is located around to = 0.14, 
which is consistent with previous findings by Sakkos et 
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FIG. 2. Influence of the relative interslice spacing to for N = 
4, r a = 4 and 9 = 0.5 on the convergence of the propagator. 
The exact result known from CPIMC (green line) is compared 
to the PB-PIMC results for P = 2, P = 3 and P = 4 for the 
fixed free parameter ai = 0.33 over the entire to range. The 
optimal value is located around to = 0.14. 


al. [45] and the application of PB-PIMC to electrons in 
a quantum dot [39]. For completeness, we mention that 
the kinetic energy I\ exhibits the same behavior. Hence, 
we use the combination ai = 0.33 and to = 0.14 for all 
presented simulations in this work. However, it should be 
noted that our method converges for all possible choices 
of the free parameters. In Fig. 3, we demonstrate the 
convergence of the energy with respect to the number of 
propagators for the same system as in Fig. 2. However, 
since V and K nearly cancel for this particular combi¬ 
nation of r s , 9 and N, we investigate the convergence 
of both contributions separately. The top panel shows 
the potential energy versus the inverse number of prop¬ 
agators P” 1 oc e and we compare the PB-PIMC results 
to the exact value (with the corresponding confidence 
interval) from CPIMC. We find that as few as two prop¬ 
agators allow for a relative accuracy AV/|V’| ~ 10~ 4 and 
with P = 4 the potential energy is converged within error 
bars. In the bottom panel, we show the same information 
for the kinetic energy K. The variance of K is one order 
of magnitude larger than that of V and, for two propaga¬ 
tors, we find the relative time step error A.K/K ~ 10 -3 . 
With increasing P, the PB-PIMC results are fluctuating 
around the exact value, within error bars. 

Finally, we adress the r s -dependence of the time step 
error by comparing PB-PIMC results for V with P = 2 
(red crosses) and P = 3 (blue squares) to the exact val¬ 
ues from CPIMC. In Fig. 4, the relative error of the po¬ 
tential energy AP/|P| is plotted versus r s for N = 4 
spin-polarized electrons at 9 = 0.5. The increased error- 
bars for larger r s are a manifestation of the sign problem 
from CPIMC [32], while for the rest the statistical un¬ 
certainty from PB-PIMC predominates. The time step 
error is smaller for three propagators over the entire r s — 
range, as it is expected, and adopts a maximum around 



1/P 

FIG. 3. Convergence of the potential (top) and kinetic (bot¬ 
tom) energy for N = 4, r s = 4 and 6 = 0.5 with to = 0.14 and 
ai = 0.33. In the top panel, the potential energy V is plot¬ 
ted versus the inverse number of propagators oc e and 
the PB-PIMC results are compared to the exact value known 
from CPIMC. The bottom panel shows the same information 
for the kinetic energy K. 


r s = 1. This can be understood by recalling the source 
of the systematic error in PB-PIMC. For r s —» 0, the 
UEG approaches an ideal system and the commutator 
error from K and V vanishes. For r s —> oo, on the other 
hand, the particles are more separated and the system 
becomes more classical. Therefore, the neglected com¬ 
mutator terms are most important at indermediate r s , 
which is the case for the results in Fig. 4. 

We conclude that as few as two or three propagators 
provide sufficient accuracy to assess the discrepancy be¬ 
tween CPIMC and RPIMC observed in previous studies 
[33]. In particular, the selected benchmark temperature, 
9 = 0.5, is even lower than for all other simulations to 
be presented in this work. Hence, the observed time step 
error constitutes an upper bound for the accuracy of our 
results in the remainder of the paper. 
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FIG. 4. Accuracy of two and three propagators over a broad 
r a range for N = 4 and 9 = 0.5 with to = 0.14 and an = 
0.33. We show the relative difference between the potential 
energy from PB-PIMC and CPIMC, AV/|V|, for the optimal 
parameters from the fourth order propagator. 



r s 

FIG. 5. The average sign is plotted versus the density pa¬ 
rameter r s for three different temperatures and N = 33 spin- 
polarized electrons with P = 2, oi = 0.33 and to = 0.14. 

B. Density parameter dependence 

Among the most interesting questions regarding the 
implementation of PB-PIMC for the UEG is the range 
of applicability with respect to the density parameter r s . 
To adress this issue, we simulate N = 33 spin-polarized 
electrons, which corresponds to a closed momentum shell 
and is often used as a starting point for finite size cor¬ 
rections. In Fig. 5, we show the average sign S versus 
r s for three different temperatures over a broad density 
range. All curves exhibit a qualitatively similar behav¬ 
ior, that is, a smooth decrease of S towards smaller r a 
until it saturates. At large r s , the coupling induced par¬ 
ticle separation mostly exceeds the extension of the sin¬ 
gle particle wavefunctions and quantum exchange effects 


do not play a dominant role. With decreasing r s , the 
UEG approaches an ideal system and the particles be¬ 
gin to overlap, which leads to sign changes in the deter¬ 
minants. However, due the blocking, the average sign, 
instead of dropping exponentially, remains finite which 
implies that, for the three depicted temperatures, PB- 
PIMC is applicable over the entire density range. This 
is in stark contrast to standard PIMC, see e.g. supple¬ 
ment of [25]. Nervertheless, with decreasing temperature 
the sign drops and the FSP makes the simulations more 
involved, cf. section IIIC. 

In Fig. 6, we compare the corresponding energies with 
RPIMC [49] and CPIMC [33], where they are available. 
The top row displays the relative difference in the po¬ 
tential energy towards PB-PIMC with two propagators. 
For 6 = 4: and 6 = 2, we find excellent agreement with 
CPIMC. For the lowest temperature, 9 = 1 , the CPIMC 
values are systematically lower by AU/|U| < 10 -3 . How¬ 
ever, this discrepancy can be explained by the conver¬ 
gence behavior of the propagator, cf. Fig. 4, since the 
potential (and kinetic) energy is expected to converge 
from above towards the exact result. To confirm this 
assumption, we also plot results for P = 3 and 9=1 , 
visualized by the grey triangles. Evidently, these points 
coincide with the CPIMC data everywhere within the er- 
rorbars and, thus, can be regarded as quasi-exact. The 
RPIMC data for V, on the other hand, exhibit a system¬ 
atic discrepancy with respect to PB-PIMC and CPIMC 
[33]. At r s = 1, the energies approximately differ by 
AU/|U| ~ 0.02, but the difference decreases with increas¬ 
ing r s . In the center row, we display the relative differ¬ 
ence in the kinetic energy. Again, all PB-PIMC results 
are in good agreement with CPIMC. On the other hand, 
there is no clear systematic deviation between the PB- 
PIMC and RPIMC data, although most RPIMC-values 
for 9 = 1 are lower while the opposite holds for most 
values for 9 = 4. Finally, the bottom row displays the 
relative difference in the total energy. Interestingly, for 
9 = 1 the difference of RPIMC in V and I\ towards PB- 
PIMC nearly cancels, so that E appears to be in good 
agreement. In particular, even the value for 9 = 1 and 
r s = 4, where the potential energy is an outlier, and 
both V and K exhibit a maximum deviation, is almost 
within single error bars. For completeness, we have also 
included the total energy for 9=1 and r s = 40 from 
standard PIMC [49], cf. the black cirlce, which is in ex¬ 
cellent agreement with PB-PIMC as well. For 9 = 2 and 
9 = 4, most RPIMC values for E are higher than PB- 
PIMC, although the deviation hardly exceeds twice the 
error bars. 


C. Temperature dependece 

Finally, we investigate the performance of PB-PIMC 
with respect to the temperature. In Fig. 7, the average 
sign is plotted versus 9 for N = 33 spin-polarized elec¬ 
trons at r s = 10, r s = 1 and r s = 0.1. All three curves 
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FIG. 6. Comparison of PB-PIMC with CPIMC and RPIMC for N = 33 spin-polarized electrons and three temperatures. In 
the top row, the relative deviation of the potential energy towards PB-PIMC with P = 2, to = 0.14 and a\ = 0.33 is plotted 
versus r s . The center and bottom rows display the same information for the kinetic and total energy, respectively. 


exhibit a similar behavior, that is, a large sign S at high 
temperature and a monotonous decay for T —> 0. How¬ 
ever, for r 8 = 10, the system is significantly less degen¬ 
erate than for both other density parameters, and even 
at 9 = 0.5, the average sign of S ss 0.056 indicates that 
the simulations are feasible. For r s = 1 and r s = 0.1, 
the decay of S is more rapid and, at low temperature, 
the simulations are more involved. In particular, half the 
Fermi temperature seems to constitute the current limit 
down to which reasonable results can be achieved for such 
r s -values (and this particle number) and, for r s = 0.1, 
the sign is zero within error bars, cf. the dashed line. Fi¬ 
nally, we note that the average signs for the two smaller 


depicted r s parameters are more similar to each other 
than to r s = 10. We characterize the temperature in 
units of the ideal Fermi temperature, which is appropri¬ 
ate for weak coupling. However, for large r a , the system 
becomes increasingly nonideal and, therefore, 9 does not 
constitute an adequate measure for the degeneracy. 

In Fig. 8, we compare the energies of the N = 33 elec¬ 
trons at r s = 1 from PB-PIMC both to RPIMC [49] and 
CPIMC. The top panel displays the relative difference 
in the potential energy versus 9. The CPIMC results 
for V are in good agreement with PB-PIMC, while the 
RPIMC data are systematically higher, by about 2%. In¬ 
terestingly, this behavior appears to be almost indepen- 
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FIG. 7. The average sign is plotted versus the temperature 8 
for r B = 10, r s = 1 and r a =0.1 and N = 33 spin-polarized 
electrons with P = 2 and the free parameters to = 0.14 and 
ai = 0.33. 




FIG. 8. Comparison with CPIMC and RPIMC as a function 
of temperature. In the top panel, the relative deviation of 
the potential energy from the PB-PIMC result is plotted ver¬ 
sus 9 for N = 33 spin-polarized electrons and r s = 1. The 
bottom panel displays the same information for the kinetic 
contribution. 


dent of the temperature. In the bottom panel, the same 
information is shown for the kinetic energy and, again, 
PB-PIMC agrees with CPIMC over the entire tempera¬ 
ture range. The large statistical uncertainty at 9 = 0.5 is 
a manifestation of the FSP in PB-PIMC, which prevents 
us from obtaining more precise kinetic energies with fea¬ 
sible computational effort. The RPIMC data for K are 
slightly lower, at low temperature, which confirms the 
trend observed by Schoof et al. [33], and seems to con¬ 
verge towards the other methods for large 6. 


IV. DISCUSSION 

In summary, we have successfully extended the permu¬ 
tation blocking path integral Monte Carlo (PB-PIMC) 
method [39] to the uniform electron gas at finite tem¬ 
perature. We have started the discussion with a brief 
introduction of our simulation scheme, which combines 
a fourth-order factorization of the density matrix with 
the application of antisymmetric imaginary time prop¬ 
agators, i.e., determinants. This allows us to combine 
permutations, which appear as individual configurations 
with positive and negative sign in standard PIMC, into a 
single configuration weight. Therefore, the average sign 
in our simulations is significantly increased. 

To assert the quality of our numerical results, we have 
investigated the optimization of the free parameters of 
our propagator and demonstrated the convergence of 
both the potential and kinetic energy with respect to 
the number of imaginary time steps. We have found that 
even for the lowest considered temperature, 9 = 0.5, as 
few as two propagators allow for a relative accuracy of 
0.1% and 0.01% in the kinetic and potential energy, re¬ 
spectively. After this preparatory work, we have shown 
results for N = 33 spin-polarized electrons, which is a 
commonly used model system as it is well suited to be 
a starting point for the extrapolation to the macroscopic 
limit (finite size corrections). Interestingly, PB-PIMC 
is feasible over the entire density range and, therefore, 
allows us to compare our results to both CPIMC and 
RPIMC data, where they are available. Our PB-PIMC 
data exhibit a very good agreement with CPIMC, for 
both the potential and kinetic energy, for all three in¬ 
vestigated temperatures. On the other hand, we observe 
deviations between PB-PIMC and RPIMC of up to 3% 
in the potential energy, which decreases towards strong 
coupling. For the kinetic energy, we find no systematic 
trend although, for 6 = 1, most of the RPIMC-values 
are smaller while, for 9 = 4, most are larger than the 
PB-PIMC results. However, for both temperatures this 
deviation hardly exceeds twice the RPIMC errorbars. 

Finally, we have investigated the applicability of PB- 
PIMC to the N = 33 spin-polarized electrons with re¬ 
spect to the temperature. With decreasing 9, exchange 
effects lead to more negative determinants in the configu¬ 
ration weights and, therefore, a smaller average sign. For 
the physically most interesting density regime, r s ~ 1, 


































simulations are feasible above 9 = 0.5 while for larger r s 
even lower temperatures are possible. A comparison of 
the energies for r s = 1 over the entire applicable temper¬ 
ature range has again revealed an excellent agreement 
with CPIMC. On the other hand, we observe a nearly 
^-independent relative deviation between PB-PIMC and 
RPIMC in the potential energy of approximately 2%, 
whereas differences in the kinetic energy are observed 
only towards low temperature. 

We conclude that our permutation blocking PIMC ap¬ 
proach is capable to provide accurate results for the UEG 
over a broad parameter range. This approach is efficient 
above a minimum temperature of about 0.5 Tp and, thus, 
complements CPIMC. Even though PB-PIMC carries a 
small systematic error (which is controllable and depends 
only on the number of time slices), we expect it to be use¬ 
ful for the development and test of other new techniques 
such as DMQMC [34, 35] and other novel versions of 
fermionic PIMC, such as the approximate treatment of 
exchange cycles by DuBois et al. [50] or a variational ap¬ 
proach to the RPIMC nodes, e.g. [51]. 

A natural follow-up of this work will be the extension of 
PB-PIMC to unpolarized systems which, in combination 
with CPIMC, should allow for a nearly complete descrip¬ 
tion of the finite temperature UEG over the entire density 


range. In addition, we aim for the application or deriva¬ 
tion of finite size corrections in order to extrapolate our 
results to the macroscopic limit [47, 52, 53] which could 
be followed by the construction of a new analytical fit 
formula for the UEG at finite temperature, e.g. [54, 55]. 
Finally, since PB-PIMC allows for efficient simulations 
in the warm dense matter regime, applications to two- 
component plasmas, such as dense hydrogen [56-58], are 
within reach. 
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